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Abstract. The statistical properties of the bid-ask spread of a frequently traded Chinese stock listed on 
the Shenzhen Stock Exchange are investigated using the limit-order book data. Three different definitions 
of spread are considered based on the time right before transactions, the time whenever the highest buying 
price or the lowest selling price changes, and a fixed time interval. The results are qualitatively similar no 
matter linear prices or logarithmic prices are used. The average spread exhibits evident intraday patterns 
consisting of a big L-shape in morning transactions and a small L-shape in the afternoon. The distributions 
of the spread with different definitions decay as power laws. The tail exponents of spreads at transaction 
level are well within the interval (2, 3) and that of average spreads are well in line with the inverse cubic 
law for different time intervals. Based on the detrended fluctuation analysis, we found the evidence of long 
memory in the bid-ask spread time series for all three definitions, even after the removal of the intraday 
pattern. Using the classical box-counting approach for multifractal analysis, we show that the time series 
of bid-ask spread does not possess multifractal nature. 

PACS. 89.65.Gh Economics; econophysics, financial markets, business and management - 89. 75. Da Sys- 
tems obeying scaling laws - 05.45.Df Fractals 
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The continuous double auction (CDA) is a dominant mar- 
ket mechanism used to store and match orders and to fa- 
cilitate trading in most modern equity markets PQ . In most 
of the order driven markets, there are two kinds of basic 
orders, called market orders and limit orders. A market or- 
der is submitted to buy or sell a number of shares at the 
market quote which results in an immediate transaction, 
while a limit order is placed to buy (or sell) a number of 
shares below (or above) a given price. All the limit orders 
that fail to result in an immediate transaction are stored 
in a queue called limit-order book. Buy limit orders are 
called bids while sell limit orders are called asks or offers. 
Best bid price b(t) and best ask (or best offer) price a(t) 
are the highest buying price and the lowest selling price 
at any time t in the limit-order book. The best bid (or 
ask) is called the same best for buy (or sell) orders, while 
the best ask (or bid) is called the opposite best for buy (or 
sell) orders. A limit order causes an immediate transac- 
tion if the associated limit price penetrates the opposite 
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best price. Such kind of limit orders are called marketable 
limit orders or effective market orders and other limit or- 
ders are termed effective limit orders. In the Chinese stock 
market, only limit orders were permitted in the placement 
of orders before July 1, 2006. 

It is a dynamic process concerning the limit-order book. 
Effective limit orders accumulate in the book while effec- 
tive market orders cause transactions and remove the limit 
orders according to their price and the time they arrive. 
Effective limit orders can also be removed by cancelation 
for a variety of reasons. Unveiling the dynamics of order 
placement and cancelation will deepen our understanding 
of the microscopic mechanism of price formation and allow 
us to reproduce remarkably many key features of common 
stocks such as the probability distribution of returns [HO 

The difference between best-ask price and best-bid price, 
s(t) = a(t) — b(t), is the bid-ask spread. Numerous work 
has been carried out to explore the different components 
of the bid-ask spread [TOllllj . On the other hand, there 
are several groups studying the statistical properties of 
the bid-ask spread time series for different stock markets. 
Farmer et al. reported that the bid-ask spread defined by 
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ln[a(i)] — ln[b(i)] on the London Stock Exchange follows 
power-law distribution in the tail 

P(> S )~ S -C, (1) 

where the exponent ( = 3.03 ± 0.41 ranging from 2.4 to 
3.9 p~2l[9] . which is well consistent with the inverse cubic 
law p~3 | fl~4j J 3] . In addition, Mike and Farmer found that 
the spread possesses long memory with the Hurst index 
being 0.75 < H < 0.85 0. Plerou et al. adopted the 116 
most frequently traded stocks on the New York Stock Ex- 
change over the two-year period 1994-1995 to investigate 
the coarse-grained bid-ask spread over a time interval At 
and found that the tail distribution decays as a power law 
with a mean tail exponent of £ = 3.0 ±0.1 and the spread 
after removing the intraday pattern exhibits long memory 
with H = 0.73 ± 0.01 [TS]. Qualitatively similar results 
were found by Cajueiro and Tabak in the Brazilian equity 
market where the mean tail exponent is £ = 2.18 ranging 
from 1.18 to 2.97 and the Hurst index is H = 0.68 ± 0.08 
varying from 0.52 to 0.89 [TB]. 

Due to the fast development of the economy of China 
and the increasing huge capitalization of its stock market, 
more concerns are attracted to study the emerging Chi- 
nese stock market. In order to reduce the market risks and 
speculation actions, the Chinese stock market adopts t + 1 
trading system, which does not allow traders to sell and 
buy stocks on the same day, and no market orders were 
permitted until July 1, 2006, which may however consume 
the liquidity of the market and cause the spread to show 
different properties when compared to other stock mar- 
kets. In this work, we investigated the probability distri- 
bution, long memory, and presence of multifractal nature 
of the bid-ask spread using limit-order book data on the 
Shenzhen Stock Exchange (SSE) in China. 

The rest of this paper is organized as follows. In Sec.[H 
wc describe in brief the trading rules of the Shenzhen 
Stock Exchange and the database we adopt. Section [3] in- 
troduces three definitions of the bid-ask spread and inves- 
tigates the intraday pattern in the spread. The cumulative 
distributions of the spreads for different definitions are dis- 
cussed in SecHl We show in Sec.[S]the long memory of the 
spread based on the detrended fluctuation analysis (DFA) 
quantified by the estimate of the Hurst index. In Sec.[6j we 
perform multifractal analysis on the bid-ask spread time 
scries. The last section concludes. 



2 SSE trading rules and the data set 

Our analysis is based on the limit-order book data of a 
liquid stock listed on the Shenzhen Stock Exchange. SSE 
was established on December 1, 1990 and has been in op- 
eration since July 3, 1991. The securities such as stocks, 
closed funds, warrants and Lofs can be traded on the Ex- 
change. The Exchange is open for trading from Monday to 
Friday except the public holidays and other dates as an- 
nounced by the China Securities Regulatory Commission. 
With respect to securities auction, opening call auction 



is held between 9:15 and 9:25 on each trading day, fol- 
lowed by continuous trading from 9:30 to 11:30 and 13:00 
to 15:00. The Exchange trading system is closed to orders 
cancelation during 9:20 to 9:25 and 14:57 to 15:00 of each 
trading day. Outside these opening hours, unexecuted or- 
ders will be removed by the system. During 9:25 to 9:30 
of each trading day, the Exchange is open to orders rout- 
ing from members, but docs not process orders or process 
cancelation of orders. 

Auction trading of securities is conducted either as a 
call auction or a continuous auction. The term "call auc- 
tion" (from 9:15 to 9:25) refers to the process of one-time 
centralized matching of buy and sell orders accepted dur- 
ing a specified period in which the single execution price is 
determined according to the following three principles: (i) 
the price that generates the greatest trading volume; (ii) 
the price that allows all the buy orders with higher bid 
price and all the sell orders with lower offer price to be 
executed; and (iii) the price that allows either buy orders 
or sell orders to have all the orders identical to such price 
to be executed. 

The term "continuous auction" (from 9:25 to 11:30 and 
from 13:00 to 15:00) refers to the process of continuous 
matching of buy and sell orders on a one-by-one basis and 
the execution price in a continuous trading is determined 
according to the following principles: (i) when the best ask 
price equals to the best bid price, the deal is concluded at 
such a price; (ii) when the buying price is higher than the 
best ask price currently available in the central order book, 
the deal is concluded at the best ask price; and (iii) when 
the selling price is lower than the best bid price currently 
available in the central order book, the deal is executed 
at the best bid price. The orders which arc not executed 
during the opening call auction automatically enter the 
continuous auction. 

The tick size of the quotation price of an order for A 
shared is RMB 0.01 and that for B sharctQ is HKD 0.01. 
Orders are matched and executed based on the principle 
of price-time priority which means priority is given to a 
higher buy order over a lower buy order and a lower sell 
order is prioritized over a higher sell order; The order se- 
quence which is arranged according to the time when the 
Exchange trading system receives the orders determines 
the priority of trading for the orders with the same prices. 

We studied the data from the limit-order book of the 
stock SZ000001 (Shenzhen Development Bank Co., LTD) 
in the whole year of 2003. The limit-order book recorded 
high-frequency data whose time stamps are accurate to 



1 A shares are common stocks issued by mainland Chinese 
companies, subscribed and traded in Chinese RMB, listed in 
mainland Chinese stock exchanges, bought and sold by Chinese 
nationals. A-share market was launched in 1990. 

2 B shares are issued by mainland Chinese companies, 
traded in foreign currencies and listed in mainland Chinese 
stock exchanges. B shares carry a face value denominated in 
Renminbi. The B Share Market was launched in 1992 and was 
restricted to foreign investors before February 19, 2001. B share 
market has been opened to Chinese investors since February 
19, 2001. 
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0.01 second. The size of the data set is 3, 925, 832, includ- 
ing 12, 965 invalid orders, 122, 034 order submissions and 
cancelations in the opening call auction, 47, 576 order sub- 
missions and cancelations during the cooling period (9:25- 
9:30), and 3,743,257 valid events during the continuous 
auction. In continuous auction, there are 317, 015 cancela- 
tions of buy orders and 274, 929 cancelations of sell orders, 
889,700 effective market orders, and 2,261,613 effective 
limit orders. Table [1] shows a segment taken from the limit- 
order book recorded on 2003/07/09. The seven columns 
stand for order size, limit price, time, best bid, best ask, 
transaction volume, and buy-sell identifier (which identi- 
fies whether a record is a buy order, sell order, or a cance- 
lation). For a cancelation record, the limit price is set to 
be zero. 



Table 1. A segment of the limit-order book 
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3 Defining bid-ask spread 

The literature concerning the bid-ask spread gives differ- 
ent definitions [6 T l9 [ ll0 p ii p i2 p i5 p i6 p i7 p i8 p i9j. In this sec- 
tion, we discuss three definitions according to sampling 
time when best bid prices and best ask prices are selected 
to define the spread. Some definitions are based on the 
transaction time, while the others are based on the phys- 
ical time. The latter scheme is actually a coarse-graining 
of the data within a given time interval. 

3.1 Definition I 

The first definition of the bid-ask spread used in this work 
is the absolute or relative difference between the best ask 
price and the best bid price right before the transaction, 
that is, 

s(t) = a(t) - b(t) (2a) 
for absolute difference or 

fl(t)=log 10 [a(t)]-log 10 [&(*)] (2b) 

for relative difference. This was used to analyze the stocks 
on the London Stock Exchange [12, 9J. The size of the 
spread time series is 895, 606. 

3.2 Definition II 

The best ask price or the best bid price may change due 
to the removal of all shares at the best price induced by an 



effective market order, or the placement of an limit order 
inside the spread, or the cancelation of all limit orders at 
the best bid/ask price. Hence the bid-ask spread does not 
always change when a transaction occurs, and it never- 
theless changes without transaction. This suggests to in- 
troduce an alternative definition of the spread which con- 
siders the absolute or relative difference between the best 
bid price and the best ask price whensoever it changes. 
The expressions of definition II are the same as those in 
Eq. ((2]) except that they have different definitions for the 
time t. The size of the spread time series is 142, 913. 

3.3 Definition III 

Obviously, the time in the first two definitions are on the 
basis of "event" . An alternative definition considers the 
average bid-ask spread over a time interval At [2D]. In 
this definition, the bid-ask spread is the average difference 
between the best ask the best bid when transactions occur 
over a fixed time interval [15] : 

1 N 

S W = ' *(*) = ~ M*) > ( 3 ) 

i=l 

where ai{t) and bi(t) are the best ask and bid prices in 
the time interval (t — At, t], and N is the total number of 
transaction in the interval and is a function of t and At. 
We use At = 1, 2, 3, 4, and 5 minute(s) to calculate the 
average spreads. 

3.4 Intraday pattern 

In most modern financial markets, the intraday pattern 
exists extensively in many financial variables [2T1I22P23] . 
including the bid-ask spread [20] . The periodic pattern 
has significance impact on the detection of long memory 
in time series [23]. To the best of our knowledge, the inves- 
tigation of the presence of intraday pattern in the spreads 
of Chinese stocks is lack. 

Figurc[T]shows the autocorrelation function (s(t)s(t + £)} 
as a function of the time lag I for the average bid-ask 
spread calculated from definition III with linear best bids 
and asks. We note that the results are very similar when 
logarithmic prices are adopted in the definition. We see 
that there are spikes evenly spaced along multiples of 245 
min, which is exactly the time span of one trading day. 
What is interesting is that Fig. Q] indicates that the aver- 
age spread also possesses half-day periodicity. 

In order to quantify the intraday pattern, we introduce 
a variable A(t), which is defined as the average bid-ask 
spread at time t for all the trading days, that is, 

M 

^) = mE s ^)> ( 4 ) 

where M is the number of trading days in the data set and 
s J (i) is the bid- ask spread at time t of day j. The spread 
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Fig. 1. Autocorrelation function {s(t)s(t + £)} of the average 
bid-ask spread calculated from definition III with the time in- 
terval At — 1 min. Note that one trading day contains 245 
trading minutes in the Chinese stock market. 



S(t) after removing the intraday pattern reads 
S(t) = s(t)/A(t) . 



(5) 



Figure [2] illustrates the intraday pattern of the bid-ask 
spread with At = 1 minute. The overall plot shows an 
evident L-shaped pattern, which is consistent with the 
one-day periodicity shown in the autocorrelation function 
in Fig. [TJ After the opening call auction, the spread A(t) 
widens rapidly and reaches its maximum 0.0183 at the end 
of the cooling auction (9:30 Then it decreases sharply 
in fifteen minutes and becomes flat at a level of 0.0112 ± 
0.0008 afterwards till 11:30. At the begin of continuous 
auction in the afternoon, A(t) abruptly rises to 0.0133 
and drops down to a stable level within about ten minutes 
which maintains until the closing time 15:00. Therefore, 
there are two L-shaped patterns each day, which suggests 
that the wide spread is closely related to the opening of the 
market. The intraday pattern makes no difference when we 
use At = 2, 3, 4, and 5 minutes. 



4 Probability distribution 

The cumulative distributions of the bid-ask spread of stocks 
in different stock markets decay as power laws with the 
tail exponent close to 3 for the major western markets 
[12, 15, 9J and much smaller and more heterogeneous in 
an emerging market |16j . Similar behavior is found in the 
Chinese stock market. Figure [3] presents the complemen- 
tary cumulative distribution s) of the spreads using 
definition I and II, where linear prices are used. Since the 



3 In 2003, there were three best prices at each side disposed 
in 9:25 and remained unchanged during the cooling period. 
Hence, the spreads shown in Figure [2] during this period are 
virtually genrated according to the trading mechanism. 
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Fig. 2. Intraday pattern in the bid-ask spread with At = 1 
min. The spread reaches its maximum at the end of the cooling 
period at 9:30. 



minimum spread equals to the tick size 0.01, the abscissa 
is no less than -2 in double logarithmic coordinates and 
0.01) = 1 for both definitions. The proportion of s = 
0.01 in the first definition is much greater than in the sec- 
ond definition such that the s) for the second defin- 
tion drops abruptly for small spreads s. The two distribu- 
tions decay as power laws with exponents (j = 2.57 ±0.06 
for definition I and (n = 2.30 ± 0.05 for definition II. 
When logarithmic prices are utilized, the spreads also fol- 
low power-law tail distributions with (j = 2.67 ± 0.03 for 
definition I and £n = 2.42±0.04 for definition II. Not much 
difference in the corresponding tail exponents (j and (n 
was found for logarithmic and linear prices. 
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Fig. 3. Empirical complementary cumulative distribution of 
the spreads calculated from definitions I and II using linear 
prices. 



Figure 3] illustrates the complementary cumulative dis- 
tributions of the average spreads over time interval At = 
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1, 2, 3, 4, and 5 minute(s) calculated from definition III 
with linear prices. The average spreads have power-law 
tails with the exponents equal to Cm,i = 2.99 ± 0.04, 
Cm,2 = 3.00±0.04, Cm,3 = 3.00±0.05, Cm,4 = 2.95±0.05, 
and Cm, 5 = 2.97 ± 0.06. Similarly, for logarithmic prices, 
we find similar power-law tail distributions with Ciu.i = 
3.07 ± 0.06, Cm 2 = 2.95 ± 0.05, Cm 3 = 3.00 ± 0.04, 
Cm,4 = 2.97 ± 0.07, and Cm, 5 = 2.98 ± 0.07. We find that 
all the tail exponents Cm, At for both linear and logarith- 
mic prices are very close to three and are independent to 
the time interval St, showing a nice inverse cubic law. This 
is well in agreement with the results in the NYSE case for 
At = 15, 30, and 60 min [15]. 




-2 -1.8 -1.6 -1.4 -1.2 -1 



log 10 s 

Fig. 4. Empirical complementary cumulative distributions of 
the average spreads calculated from definition III with time 
intervals At = 1, 2, 3, 4, and 5 min using linear prices. The 
markers represent the real data and the solid lines are the best 
fits in the scaling ranges. The curves with At > 1 has been 
translated vertically for clarity. 



There arc also significant discrepancies. Comparing the 
cumulative distributions in Fig. 2] and that on the NYSE 
[15] , significant differences are observed. The distribution 
of the spreads on the SSE decays much faster than that 
on the NYSE for small spreads. In other words, the pro- 
portion of small spreads is much larger on China's SSE. 
Possible causes include the absence of market orders, no 
short positions, the maximum percentage of fluctuation 
(10%) in each day, and the t+ 1 trading mechanism in the 
Chinese stock markets on the one hand and the hybrid 
trading system containing both specialists and limit-order 
traders in the NYSE on the other hand. The exact cause is 
not clear for the time being, which can however be tested 
when new data are available after the introduction of mar- 
ket orders in July 1, 2006. Moreover, the PDF's in SSE 
drop abruptly after the power-law parts for the largest 
spreads, which is not observed in the NYSE case [15] . 



5 Long memory 

Another important issue about financial time series is the 
presence of long memory, which can be characterized by 
its Hurst index H. If H is significantly larger than 0.5 the 
time series is viewed to possess long memory. Long mem- 
ory can be defined cquivalently through autocorrelation 
function C{£) ~ £~ 7 and the power spectrum p(uj) ~ w~ T ', 
where the autocorrelation exponent 7 is related to the 
Hurst index H by 7 = 2 — 2H [25,26], and the power 
spectrum exponent ?/ is given by r) = 2H — 1 [2"Tll28| . 

There are many methods proposed for estimating the 
Hurst index such as the rescaled range analysis (RSA) [29] 
l30ll3Tl[32l[33ll34] . fluctuation analysis (FA) [35], detrended 
fluctuation analysis (DFA) [36ll24lf25] . wavelet transform 
module maxima (WTMM) method [37135113211101111], de- 
trended moving average (DMA) [42ll43U44l[45lf46] . to list a 
few. We adopt the detrended fluctuation analysis. 

The method of detrended fluctuation analysis is widely 
used for its easy implementation and robust estimation 
even for a short time series [47 ] l48 p 28 p 49] . The idea of 
DFA was invented originally to investigate the long-range 
dependence in coding and noncoding DNA nucleotides 
sequence [36] and then applied to various fields including fi- 
nance. The method of DFA consists of the following steps. 

Step 1: Consider a time series x(t), t — 1, 2, ■ • • , N. We 
first construct the cumulative sum 

t 

u(t)=J>(<), t = l,2,---,N. (6) 

Step 2: Divide the series u(t) into Ng disjoint segments 
with the same length I, where Ne — [N/l]. Each segment 
can be denoted such that u v (i) = u{l + i) for 1 ^ 

i ^ £, and I = (v — 1)1. The trend of u v in each segment 
can be determined by fitting it with a linear polynomial 
function u v . Quadratic, cubic or higher order polynomials 
can also be used in the fitting procedure while the simplest 
function could be linear. In this work, we adopted the 
linear polynomial function to represent the trend in each 
segment with the form: 

u v (i) = ai + b , (7) 

where a and b are free parameters to be determined by 
the least squares fitting method and 1 ^ i ^ I. 

Step 3: We can then obtain the residual matrix e v in 
each segment through: 

e v (i) = u v (i) - u v (i) , (8) 

where 1 ^ i ^ I. The detrended fluctuation function 
F(v,£) of the each segment is defined via the sample vari- 
ance of the residual matrix e v as follows: 

^MH^EM*)] 2 . (9) 
i=i 

Note that the mean of the residual is zero due to the de- 
trending procedure. 
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Step 4: Calculate the overall dctrcndcd fluctuation func- 
tion F(£), that is, 

V — 1 

Step 5: Varying the value of £, we can determine the 
scaling relation between the detrended fluctuation func- 
tion F{£) and the size scale £, which reads 

F(£)~l H , (11) 

where H is the Hurst index of the time series [47l[25] . 

Figure [5] plots the detrended fluctuation function F(£) 
of the bid-ask spreads from different definitions using lin- 
ear prices. The bottom F(£) curve is for the average spread 
after removing the intraday pattern. All the curves show 
evident power-law scaling with the Hurst indexes H\ = 
0.91 ± 0.01 for definition I, H u = 0.92 ± 0.01 for defini- 
tion II, Fin = 0.75 ± 0.01 for definition III, and H m = 
0.77 ± 0.01 for definition without intraday pattern, re- 
spectively. Quite similar results are obtain for logarith- 
mic prices where H\ = 0.89 ± 0.01 for definition I, Hu = 
0.91 ± 0.01 for definition II, H m = 0.77 ± 0.01 for defini- 
tion III, and Hm = 0.76 ± 0.01 for definition III without 
intraday pattern. The two Hurst indexes for definitions I 
and II are higher than their counterparts on the London 
Stock Exchange where "even time" is adopted [9]. It is 
interesting to note that the presence of intraday pattern 
does not introduce distinguishable difference in the Hurst 
index and the two indexes for definition III are also very 
close to those of average spreads in the Brazilian stock 
market and on the New York Stock Exchange where real 
time is used [TMTB] . Due to the large number of data used 
in the analysis, we argue that the bid-ask spreads investi- 
gated exhibit significant long memory. 




Fig. 5. Detrended fluctuation function F(£) for the spreads 
obtained from three definition with linear prices. The curves 
have been shifted vertically for clarity. 



6 Multifractal analysis 

In this section, we investigate whether the time series of 
bid-ask spread obtained from definition III possesses mul- 
tifractal nature. The classical box-counting algorithm for 
multifractal analysis is utilized and described below [50] . 

Consider the spread time series S(t), t = 1, 2, ■ • ■ ,N. 
First, we divide the series S(t) into Ne disjoint segments 
with the same length £, where Ng = [N /£]. Each segment 
can be denoted as S v such that S v (i) = S(l + i) for 1 ^ 
i ^ £, and I = (v — 1)£. The sum of S v in each segment is 
calculated as follows, 

i 

r(v,£) = Y / Sv(i), v = l,2,---,N*. (12) 

We can then calculate the gth order partition function 
r(q; £) as follows, 

N e 

r(q;£)=Y^[r(v,£W . (13) 

Varying the value of £, we can determine the scaling rela- 
tion between the partition function r(q; £) and the time 
scale £, which reads 

r(q;£) ~ . (14) 

Figure [6] illustrates the power-law scaling dependence 
of the partition function r(q; £) of the bid-ask spreads af- 
ter removing the intraday pattern in definition III for dif- 
ferent values of q, where both linear prices and logarithmic 
prices are investigated. The continuous lines are the best 
linear fits to the data sets. The collapse of the data points 
on the linear lines indicates evident power-law scaling be- 
tween r(q; £) and £. The slopes r(q) of the fitted lines 
are r(-4) = -5.02 ± 0.01, r(-2) = -3.01 ± 0.01, r(0) = 
-1.01 ± 0.01, r(2) = 0.99 ± 0.01, and r(4) = 2.98 ± 0.01 
for logarithmic prices and r(— 4) = —5.02 ±0.01, r(— 2) = 
-3.01 ± 0.01, r(0) = -1.01 ± 0.01, t(2) = 0.99 ± 0.01, 
and r(4) = 2.98 ± 0.01 for linear prices. We notice a nice 
relation r(q) = q — 1. 

Quantitatively similar results are obtained when the 
intraday pattern is not removed. The scaling exponents 
are r(-4) = -5.03 ± 0.01, r(-2) = -3.01 ± 0.01, r(0) = 
-1.01 ± 0.01, t(2) = 0.99 ± 0.01, and r(4) = 2.96 ± 0.01 
for logarithmic prices and r(— 4) = —5.03 ±0.01, r(— 2) = 
-3.02±0.01,t(0) = -1.01 ± 0.01, r(2) = 0.99 ±0.01, and 
r(4) = 2.97±0.01 for linear prices. Again, we observe that 
r(q) =q-l. 

In the standard multifractal formalism based on parti- 
tion function, the multifractal nature is characterized by 
the scaling exponents r(q). It is easy to obtain the general- 
ized dimensions D q — r(q)/(q — 1) 51,52,53 and the sin- 
gularity strength function a(q), the multifractal spectrum 
/(a) via the Legendre transform [50] : a(q) = dr(g)/dg 
and f(q) =qa- r(q). 

Figure [7] shows the multifractal spectrum /(a) and the 
scaling function r(q) in the inset for linear and logarith- 
mic prices. One finds that the two r(q) curves are linear 
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Fig. 6. Log-log plots of the partition function r(q;£) of the 
bid-ask spreads calculated from definition III with the intraday 
pattern removed for five different values of q. Both linear and 
logarithmic prices are investigated (shown in the legend). The 
markers stand for the results calculated from the real data and 
the continuous lines are the best fits. The plots for q — —2, 0, 2, 
and 4 are shifted upwards for clarity. 



and r{q) = q — 1, which is the hallmark of the presence 
of monofractality, not multifractality. The strength of the 
multifractality can be characterized by the span of sin- 
gularity Aa = a max — a m in- If Aa is close to zero, the 
measure is almost monofractal. The maximum and min- 
imum of a can be reached when q — > ±oo, which can 
not be achieved in real applications. However, Aa can 
be approximated with great precision with mediate values 
of q. The small value of Aa shown in Fig. [7] indicates a 
very narrow spectrum of singularity. Indeed, one sees that 
f(a) ~ 1 and a pa 1 for all values of q. We thus conclude 
that there is no multifractal nature in the bid-ask spread 
investigated. 



7 Conclusions 

The bid-ask spread defined by the difference of the best 
ask price and the best bid price is considered as the bench- 
mark of the transaction cost and a measure of the mar- 
ket liquidity. In this paper, we have carried out empirical 
investigations on the statistical properties of the bid-ask 
spread using the limit-order book data of a stock SZ000001 
(Shenzhen Development Bank Co., LTD) traded on the 
Shenzhen Stock Exchange within the whole year of 2003. 
Three different definitions of spread are considered based 
on event time at transaction level and on fixed interval of 
real time. 

The distributions of spreads at transaction level decay 
as power laws with tail exponents well below 3. In contrast 
the average spread in real time fulfils the inverse cubic law 
for different time intervals At = 1, 2, 3, 4, and 5 min. We 
have performed the detrended fluctuation analysis on the 
spread and found that the spread time series exhibits evi- 
dent long-memory, which is in agreement with other stock 
markets. However, an analysis using the classic textbook 
box-counting algorithm does not provide evidence for the 
presence of multifractality in the spread time series. To 
the best of our knowledge, this is the first time to check 
the presence of multifractality in the spread. 

Our analysis raises an intriguing open question that is 
not fully addressed. We have found that the spread pos- 
sesses a well-established intraday pattern composed by a 
large L-shape and a small L-shape separated by the noon 
closing of the Chinese stock market. This feature will help 
to understand the cause of the wide spread at the opening 
of the market, which deserves further investigation. 

We are grateful to Dr. Tao Wu for fruitful suggestions. This 
work was partially supported by the National Natural Science 
Foundation of China (Grant No. 70501011) and the Fok Ying 
Tong Education Foundation (Grant No. 101086). 
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Fig. 7. Multifractal function f(a) of the spreads in definition 
III with the intraday pattern removed. Inset: Scaling exponents 
r(q) of partition functions as a function of q. For clarity, the 
r(q) curve for logarithmic price is shifted upwards by 1. 
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